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Nitrogen productivity, expressed as annual yield of foliage per unit of nitrogen in the foliage, is recommended for 
interpreting responses to changes in nitrogen availability. This index can be directly coupled to growth and is well suited for 
application in mathematical analyses of tree or forest growth dynamics. This paper demonstrates that the nitrogen productivity 
decreases predictably as biomass increases probably as a result of the reduction in trradiance as the canopy becomes more 
dense. For five species of conifers from northern latitudes (Picea abies Karst., Pseudotsuga menziesii Mirb, Franco, Pinus 
nigra var. maritima (Ait.} Melv., Pinus resinosa Ait., and Pinus sylvestris L.) nitrogen productivity is shown to decrease 
linearly with increasing needle biomass. For a particular species, these regressions allow for estimation of maximum needle 
biomass. 
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Pour interpréter les réactions au changement de l'azote disponible, on propose la notion de productivité de l'azote, exprimée 
comme une production annuelle de feuillage par unité d’azote dans le feuillage. Cet indice peut être associé directement à la 
production et s'applique très bien a des analyses mathématiques de dynamique d’accroissement d’arbres ou de forêts, L'auteur 
montre que la productivité de l'azote décroit quand la biomasse s'accroît, probablement à cause de la réduction de rayonnement 
avec un couvert de plus en plus dense. Chez cing essences résineuses des régions nordiques (Picea abies Karst., Pseudotsuga 
menziesii (Mirb.) Franco, Pinus nigra var. maritima (Ait.) Melv., Pinus resinosa Ait. et Pinus sylvestris L.), la productivité 
de l'azote décroit de façon linéaire avec l'accroissement de la biomasse du feuillage. Pour une essence donnée, ces régressions 


permettent d’évaluer la biomasse maximale du feuillage. 


Introduction 


Traditionally, studies of tree nutrition have focused 
upon changes in concentration of various nutrients in 
different tissues (cf. Morrison 1974). Nutrient con- 
centrations, particularly of the foliage, have served as 
an index of the nutritional status of trees or stands. 
However, observed variations in nutrient concentra- 
tions have not been easily associated with produc- 
tion. One difficulty is that the same amount of nutrients 
can be allocated into variable amounts of foliage so 
that it might be better to express results as the product 
of concentration times tissue. This represents the 
“amount” of nutrient committed to growth of foliage. 
In this paper, I construct a framework for analysis, 
called “nitrogen productivity,” based on this principle 
and relate it to growth. 

I treat only nitrogen because it is commonly one of 
the most growth-limiting macronutrients and data are 
therefore readily available. I propose that the amount of 
biomass produced is directly related to the amount of 
nitrogen in that biomass through the nitrogen pro- 
ductivity. For field applications nitrogen productivity is 
best described as a collective property of a stand rather 
than for single trees. For simplicity, I restrict the anal- 
ysis to the production of foliage, recognizing that the 
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growth of other tissue can be estimated from correla- 
tions with this variable (e.g., Ovington 1957; Whittaker 
and Woodwell 1968; Tamm 1974), 

The approach is designed to provide a tool for 
analyses of rapidly changing growth of forests as 
well as other ecosystems. I therefore stress this objec- 
tive. Other emphases might warrant, of course, other 
analyses. 


Theory 


Quite generally, the net increment in foliage may be 
described as the difference between production of new 
foliage and loss of old foliage. If the production rate is 
assumed proportional to the amount of nitrogen in the 
foliage (N) and the loss rate to the weight of foliage 
(W), the net increment in foliage can be written 


[1] dW/dt = Pm N- fW 


Here, Piw is the nitrogen productivity and fis the rate 
of tissue mortality. Both depend in principle not only 
upon biomass (W ) and nitrogen (N) but also upon exter- 
nal variables such as weather and time (t). Equation 1 
defines the nitrogen productivity as the amount of bio- 
mass produced per unit of nitrogen in that biomass per 
unit of time. 

If we assume the rate of mortality (f) to be approxi- 
mately constant, then, a constant nitrogen productivity 
(Pi) would allow the foliage to increase indefinitely 
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with added amounts of nitrogen, an obvious impossibil- 
ity. Hence, the nitrogen productivity should decrease as 
foliage biomass (W ) increases. A decrease in the nitro- 
gen productivity could result from internal shading, 
water stress, or increased respiration by nonphoto- 
synthesizing tissue. Let us for simplicity assume that 
Py, will decrease linearly. Hence, 


[2] Pu =a-bWw 


where a and b are species-specific parameters. Equa- 
tions | and 2 imply that there is an upper limit to the 
amount of foliage, Wmax: 


(él f/(N/W max 


[3] W ma = b 


where (N/W )max is the maximum optimal nitrogen con- 
centration obtained without tissue injury (cf., Ingestad 
et al. 1981, Fig. 1). Equation 3 can be used for esti- 
mating the maximal leaf biomass of a forest once we 
have estimated the constituting parameters. 

How much might self-shading account for the pre- 
dicted reduction in the nitrogen productivity with in- 
creasing foliage biomass? In the Appendix, I present 
calculations demonstrating that simple application of 
Lambert—Beer’s law to light penetration through foli- 
age can fully account for the phenomenon. Derived 
from this exercise is an equation for canopy photo- 
synthesis such that: 


bo Y 
Aa asi 


where y is the light extinction coefficient (the frac- 
tion of light absorbed per unit of foliage penetrated), 
lo the light intensity above the canopy, and a the 
light intensity at which the photosynthetic rate is half 
that at saturation. 


Results 


Let us now apply these formulae to independent sets 
of data. I have chosen to select examples dealing with 
just foliage production to minimize confusion and cal- 
culations. The examples could easily be expanded to 
predict forest biomass accumulation. The foliage is well 
suited as a test tissue because, unlike the stem, it con- 
sists entirely of living tissue. 

The data base available for evaluating the hypothesis, 
Eq. 2 only permits one comparison during a season, 
although the possibility exists of applying the equation 
to intervals over the growing season. In conifers with a 
relatively stable needle biomass it is possible to define 
an annual nitrogen productivity without concern for 
large within-year variations. For convenience, | define 
the nitrogen productivity from values at the end of the 
growing season: Py, = (current needle biomass}/(total 
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nitrogen in all foliage). With deciduous species or with 
conifers with short-lived needles (e.g., Pinus taeda) 
this simple approach is not possible, and repeated 
samplings over the whole growing season would be 
required (Ingestad 1979a, 1980, 1981; Ericsson 198 1a, 
1981b). 

Adequate data for testing the relationship formu- 
lated in Eq. 2 were available for Picea abies Karst. 
(Norway spruce), Pseudotsuga menziesii Mirb. Franco 
(Douglas-fir), and three species of Pinus (P. nigra var. 
maritima (Ait.) Melv. (Corsican pine), P. resinosa Ait. 
(Red pine), and P. sylvestris L. (Scots pine)) covering 
a wide range of ages, stand densities, soil conditions, 
and geographical locations (Table |). Figure | displays 
graphically the relations, and statistical analyses are 
presented in Table 2. For Picea abies and Pinus 
sylvestris only are the data bases adequate (more than 
20 samples) for statistical analysis. The other three 
species, however, appear to exhibit similar relation- 
ships. Although many other studies have been reported 
in the literature, most lack information on the necessary 
data (current needie biomass, total needle biomass, and 
the nitrogen pool in foliage). 

Table 3 presents the results of the analysis of the 
effect of internal shading upon the nitrogen pro- 
ductivity. Jarvis et al. (1976) report values of 
0.35—0.80 ha - ha ' for the light extinction coefficient 
(y) for several coniferous species and under different 
light conditions. Data in Albrektson and Aronsson 
(1983), Miller and Miller (1976), and Ovington (1957) 
indicate that | kg of Scots pine needles have a projected 
area of around 3+ 107? ha, whereas Nihlgård (1972) 
gives 6.4 + 107 kg dry weight ha’! for Norway spruce, 
thus giving values of y between 10™* and 5+ 10™* kg 
dry weight ha~'. For further analysis y was assumed 
to be 3:10 kg dry weight ha™'. An important 
question is which value should be chosen for /ọ, the 
light intensity above the canopy. In Table 3, fy = 
400 pE - m7’ - s7! corresponding to an average value 
over the growing season is used (Linder and Lohammar 
1981). With these values, good agreement between 
the two sides of Eq. 4 is obtained. From this analysis I 
therefore draw the conclusion that the main factor 
responsible for the decrease in nitrogen productivity 
with increasing needle biomass is the increasing 
internal shading within the canopy. 

From the same data as used for Fig, | it is also 
possible to get an estimate of the mortality rate of the 
foliage by assuming it equal to the ratio of current to 
total needle biomass. This is correct if the needle bio- 
mass is in equilibrium with the environment (ceiling 
LAD but results in an overestimate of the mortality rate 
if canopy biomass or leaf area index are still increasing. 
Actually, surprisingly little variation in the rate of 
needle mortality in relation to total needle biomass was 


TABLE 1. Data about stands used in the analysis 


Species 


Picea abies 
Picea abies 
Picea abies 
Picea abies 
Picea abies 
Picea abies 


Pseudotsuga menziesii 
Pseudotsuga menziesii 
Pseudotsuga menziesii 


Pinus nigra 


Pinus resinosa 

Pinus resinosa 

Pinus sylvestris 
Pinus sylvestris 
Pinus sylvestris 
Pinus sylvestris 
Pinus sylvestris 
Pinus sylvestris 


Number of 
stands* 


Ages 
(years) 


28, 45, 47 


120 
64 
20 


Stand densities 
(trees/ha) 


2900—3900 
1625-3025 
880 
356 
1332 
nr 
nr 
nr 
16000 
2110 


2103—2638 
3460 
1131—1280 

760—5640 


2911, 1420, 845 


393 
815 
1095 


*Can also be number of different plots or treatments within an experiment. 


+The code refers to the symbols used in Fig. 1. 


Not reported. 


Location 


Reference 


Southern Sweden 
Central Sweden 
Southern Sweden 
Arkhangelsk 
Czechoslovakia 
Solling 
Washington State 
Washington State 
British Columbia 
England 


New York State 
New York State 
Central Sweden 
England 
Southern Finland 
Central Sweden 
Scotland 

Central Sweden 


Tamm (1968, 1975) 

Tamm (1974) 

Nihlgard (1972) 

Marchenko and Karlov (1962) 
Klimo (1980) 

Cole and Rapp (1981) 

Cole and Rapp (1981) 

Gessel et al. (1973) 

Webber (1977) 

Miller and Miller (1976) 
Miller et al. (1976) 

Wittwer et al. (1975) 

Fornes et al. (1970) 
Albrektson and Aronsson (1983) 
Ovington (1957, 1959) 
Mälkönen (1974) 

Bringmark (1977) 

Wright and Will (1958) 

Own data 
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Fic. 1. Representations of the relations between nitrogen 
productivity and needle biomass. Different data sets are rep- 
resented with different symbols, references are given in 
Table 1. The regressions for equations in Table 2 are also 
given. These regression lines are combined for comparison in 
the lower right-hand subpicture. The length of the regression 
lines there indicate the variability in the needle biomasses in 
the underlying data sets. 


observed, but to get a definite value I chose the upper 
end of the regression of mortality rate on total needle 
biomass. 

Ingestad (19795) found that the optimal nitrogen 
concentration of Picea abies and Pinus sylvestris 
both lie close to 2% of dry weight. Assuming this 
to be true for the other coniferous species (although 
van den Burg (1971) gives higher values for Pinus nigra 
and Pseudotsuga menziesii) the maximal needle bio- 
masses were calculated from Eq. 3 and are presented in 
Table 4. These values correspond quite well with the 
highest values of needle biomass reported in the 
literature. 


Discussion 


Although the data base was small, results indicate a 
general pattern. The three pine species from the north- 
ern temperate zone all exhibited similar regression 
equations (Fig. 1). The spruce, on the other hand, was 
different, having much lower nitrogen productivity 
(P) but also showing very little reduction in nitrogen 
productivity with increasing needle biomass. The latter 
should be expected because spruce is a shade-tolerant 
species and grows satisfactorily under low light. 
The Douglas-fir falls between these two. The two 
parameters a and b are strongly correlated (p = 0.952, 
P < 5%), suggesting that coniferous species may 
have evolved either to be highly efficient at high light 
intensities or to be moderately efficient at a wide range 
of light intensities. 

Precisions of the regressions presented, Table 2, are 


TABLE 2. Regressions of nitrogen productivity (a — b W) on needle biomass (W ) for different species 


Parameters with 95% confidence intervals 


Species a* bt Syn r? n 
Picea abies 18.443.8 0.000 377 + 0.000 288 2.92 0.251 24 
Pseudotsuga menziesii 34.4+22.1 0.001 17 + 0.002 24 4.55 0.480 5 
Pinus nigra 50.2+8.8 0.001 84 + 0.000 75 1.38 0.953 5 
Pinus resinosa 55.7£14.1 0.002 41 + 0.001 30 3.64 0.869 6 
Pinus sylvestris 41.424.5 0.002 04 + 0,000 74 4.98 0.465 38 


*Kilograms of dry weight per kilogram of nitrogen per year. 


‘Hectares per kilogram of nitrogen per year. 


TABLE 3. Numerical evaluation of Eq. 4 


y 

Species apEsm*+s' Wo/a +1) 
Picea abies 98 29.5x 107° 
Pseudotsuga menziesii 258 58.8x 10% 
Pinus nigra 178 46.2x 107° 
Pinus resinosa 189 48.1x 10%“ 
Pinus sylvestris 330 67.8x 10% 


*Calculated from Table 2. 


b/a* Reference for « 
20.5x 107° Fuchs et al. (1977) 
34.0x 107° Kiinstle (1971) 
36.7x 107° van den Driessche and Wareing (1966) 
43,3«K 107° Bourdeau and Laverick (1958) 
49.3x 107° Troeng and Linder (1982) 
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TABLE 4. Maximal needle biomasses 


Maximal needle biomass 


Loss rate 


(kg dry weight + ha™') 


Minimal nitrogen 


Species {year ') From Eq. 3 Observed concentration” 
Picea abies 0.15 28 300 20 600 0.82 
Pseudotsuga menziesii 0.17 22 000 14 100 0.49 
Pinus nigra 0.43 15 500 14 800 0.86 
Pinus resinosa 0.23 18 400 14 100 0.41 
Pinus sylvestris 0.39 10 700 10 500 0.94 


*Percentage of dry weight. This is the lowest nitrogen concentration that allows growth according to Eq. I 


not remarkable, but it must be remembered that the 
data sets represent stands from a variety of sites. 
Looking in some detail at the data, Fig. 1, one sus- 
pects that the nitrogen productivity of Scots pine in 
Scandinavia (data sets A, C, D, and F) is lower than 
in Britain (data sets B and E), which could be a result 
of longer periods of photosynthetic activity in 
Britain. However, the differences are not statistically 
significant. 

The data set A for Norway spruce and Scots pine are 
taken from fertilizer trials. In these data sets there are 
three points (two for the spruce and one for the pine) 
falling well below the regression lines. These points 
correspond to the highest level of fertilizer applications, 
applications that probably were super-optimal. On the 
other hand, it is possible that improved availability of 
water might increase the nitrogen productivity as indi- 
cated by the position of the irrigated plots from data 
set E for Scots pine. 

In this paper | have assumed that the nitrogen pro- 
ductivity can be treated as independent of nitrogen con- 
centration. From laboratory studies (Ingestad 1979a, 
1980, 1981) we know that this is not so. At low nitro- 
gen concentrations the nitrogen productivity decreases. 
However, as nitrogen concentrations increase, the 
nitrogen productivity increases rapidly up to a constant 
level. Thus, over wide ranges of nitrogen concen- 
trations the nitrogen productivity may, indeed, be 
regarded as constant. In situations where the nitrogen 
concentrations fall to very low levels, as might be the 
case in other ecosystems like grasslands, reformula- 
tions are necessary, not in the idea as such but in its 
interpretation, Eq. 2. 

The formulation of the nitrogen productivity (Pv) 
can be given a form that is easily interpretable biologi- 
cally and also mathematically tractable, Eq. 2. There- 
fore, the simple Eq. 1 becomes a well-formulated 
model, providing us with a convenient tool for ana- 
lysing tree and forest growth both under natural condi- 
tions and under influence of perturbations (Ingestad 
etal. 1981; Agren 1983a, 19830). Simplicity in formu- 
lations and interpretations have been the major criteria 
in this paper. | want to stress this objective because 


other objectives could have implied other interpreta- 
tions of the basic data (cf. Agren 1981). In particular, 
there is often a conflict between accuracy of quantita- 
tive prediction and qualitative understanding, where 
the former may require the inclusion of several 
details obscuring the transparency of the problem to 
the detriment of the latter. 

Leaf and Madgwick (1960) suggested that to “prop- 
erly interpret analytical data it is necessary to convert 
results from concentrations to contents on an areal 
basis.” The nitrogen productivity provides a tool for 
doing this and warrants further testing in predicting 
forest growth. The approach also appears promising for 
comparing species and interpreting their ecological dif- 
ferences. Studies on other nutrients and environmental 
factors may also justify testing within the general 
framework presented for assessing nitrogen constraints 
upon forest growth. 
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Appendix 
Derivation of Eq. 4 
Assume that the light response of the photosynthetic rate 
can be described by 


i 


[A] PU) = P,, Tio 


for the whole canopy, where / is light intensity, Pm is the 
photosynthetic rate per unit needle biomass at light saturation 


and a a parameter, Let p(z) be the needle biomass density at 
depth z in the canopy (W = | dz p(z)), and let further the light 
o 


extinction in the os be described by Lambert — Beer’s law 


my [rae pco) 


where Jo is the light intensity above the canopy and y a 
parameter. The photosynthetic rate of the whole canopy, 
Piol W ), can now be calculated 


RA] IZ) = h exp | 


ly + a 
lh exp {—yW} + a 


L 
= 7 d) = Pin 
[3A] Piss Í dz Pin IG) Fa In 
Equation 3A is an equation for the instantaneous rate of photo- 
synthesis in the whole canopy. In order to connect Pa with 
the nitrogen productivity, /o and P,, will be treated as con- 
stants. Then, Pia is also equal to the average photosynthetic 


rate. It is reasonable to expect that as photosynthetic rate goes 
up at constant nitrogen content, the ratio between the amounts 
of photosynthate and nitrogen increases. This should lead to 
an increase in nitrogen productivity. I express this relation- 
ship by assuming that the nitrogen productivity is propor- 
tional to the photosynthetic rate per unit needle biomass, 
Le., PolW) =k Pww(W), where k is the proportionality 
constant. Dividing both sides with their respective value at 
W = 0 removes the proportionality constant, yielding 
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It is clear that Eq. 4A is not an exact relationship because the 
right-hand side is not a linear function of W. However, for the 
present purpose and within the range of W that is relevant in 
this context, the right-hand side of Eq. 4A is approximated as 
a linear function. From the construction of Eq. 4A it is clear 
that the two sides are equal at W = 0. I estimate the slope of 
the right-hand side with its slope at W = 0, which gives 


b y. 
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This is the relationship between b/a and the photosyn- 
thesis related parameters Jo, a, and y that should hold if the 
decrease in nitrogen productivity could entirely be attributed 
to internal shading. 


